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The behavior near the initial state of the anisotropy parameter of the arbitrary type, homogeneous 
and anisotropic Bianchi models is considered in the framework of the brane world cosmological 
models. The matter content on the brane is assumed to be an isotropic perfect cosmological fluid 
obeying a barotropic equation of state. To obtain the value of the anisotropy parameter at an 
arbitrary moment an evolution equation is derived, describing the dynamics of the anisotropy as a 
function of the volume scale factor of the Universe. The general solution of this equation can be 
obtained in an exact analytical form for the Bianchi I and V types and in a closed form for all other 
homogeneous and anisotropic geometries. The study of the values of the anisotropy in the limit of 
small times shows that for all Bianchi type space-times flUed with a non-zero pressure cosmological 
fluid, obeying a linear barotropic equation of state, the initial singular state on the brane is isotropic. 
This result is obtained by assuming that in the limit of small times the asymptotic behavior of the 
scale factors is of Kasner-type. For brane worlds filled with dust, the initial values of the anisotropy 
coincide in both brane world and standard four-dimensional general relativistic cosmologies. 

PACS numbers; 04.20.Jb, 04.65. +e, 98.80.-k 



I. INTRODUCTION 

The idea that our four-dimensional Universe might be a three-brane, embedded in a higher dimensional space- 
time, has attracted much attention. According to the brane- world scenario, the physical fields in our four-dimensional 
space-time, which are assumed to arise as fluctuations of branes in string theories, are confined to the three brane. Only 
gravity can freely propagate in the bulk space-time, with the gravitational self-couplings not significantly modified. 
This model originated from the study of a single 3-brane embedded in five dimensions, with the 5D metric given 
by ds^ = e~^''y^ri^i,dx^dx'^ + dy'^ , which, due to the appearance of the warp factor, could produce a large hierarchy 
between the scale of particle physics and gravity. Even if the fifth dimension is uncompactified, standard AD gravity is 
reproduced on the brane. Hence this model allows the presence of large, or even infinite non-compact extra dimensions. 
Our brane is identified to a domain wall in a 5-dimensional anti-de Sitter space-time. 

The Randall-Sundrum model was inspired by superstring theory. The ten-dimensional Eg, x heterotic string 
theory, which contains the standard model of elementary particle, could be a promising candidate for the description 
of the real Universe. This theory is connected with an eleven-dimensional theory, M-theory, compactified on the 
orbifold R^^ x /Z2 0. In this model we have two separated ten-dimensional manifolds. For a review of dynamics 
and geometry of brane Universes see . 

The static Randall-Sundrum solution has been extended to time-dependent solutions and their cosmological prop- 
erties have been extensively studied H-|33|. In one of the first cosmological applications of this scenario, it was 
pointed out that a model with a non-compact fifth dimension is potentially viable, while the scenario which might 
solve the hierarchy problem predicts a contracting Universe, leading to a variety of cosmological problems ^12] . By 
adding cosmological constants to the brane and bulk, the problem of the correct behavior of the Hubble parameter 
on the brane has been solved by Cline, Grojean and Servant jl3i |. As a result one also obtains normal expansion 
during nucleosynthesis, but faster than normal expansion in the very early Universe. The creation of a spherically 
symmetric brane-world in AdS bulk has been considered, from a quantum cosmological point of view, with the use of 
the Wheeler-de Witt equation, by Anchordoqui, Nunez and Olsen [T^ . 

The effective gravitational field equations on the brane world, in which all the matter forces except gravity are 
confined on the 3-brane in a 5-dimensional space-time with Z2-symmetry have been obtained, by using a geometric 
approach, by Shiromizu, Maeda and Sasaki [l^ [l^ . The correct signature for gravity is provided by the brane with 
positive tension. If the bulk space-time is exactly anti-de Sitter, generically the matter on the brane is required to be 
spatially homogeneous. The contraction of the 5-dimensional Weyl tensor with the normal to the brane Ejj gives the 
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leading order corrections to the conventional Einstein equations on the brane. The four-dimensional field equations 
for the induced metric and scalar field on the world-volume of a 3-brane in the five-dimensional bulk, with Einstein 
gravity plus a self- interacting scalar field, have been derived by Maeda and Wands [l^ . 

Realistic brane-world cosmological models require the consideration of more general matter sources to describe the 
evolution and dynamics of the very early Universe. The influence of the bulk viscosity of the matter on the brane has 
been considered, for an isotropic flat Friedmann-Robertson- Walker (FRW) geometry, in 20] and, for a Bianchi type 
I geometry, in |23| . The first order rotational perturbations of isotropic FRW cosmological models have been studied 
in in and |2|. 

The general solution of the field equations for an anisotropic brane with Bianchi type I and V geometry, with perfect 
fluid and scalar fields as matter sources, has been obtained in Expanding Bianchi type I and V brane- worlds 

always isotropize. Anisotropic Bianchi type I brane-worlds with a pure magnetic field and a perfect fluid have also 
been analyzed . Limits on the initial anisotropy induced by the 5-dimensional Kaluza-Klein graviton stresses by 
using the CMB anisotropics have been obtained by Barrow and Maartens j^j. The dynamics of a flat, isotropic brane 
Universe with two-component matter source: a perfect fluid with a linear barotropic equation of state and a scalar 
field with a power-law potential has been investigated in 34] . Solutions for which the scalar field energy density scales 
as a power-law of the scale factor (so called scaling solutions) have been obtained and their stability analysis provided. 

A family of Bianchi type braneworlds with anisotropy has been constructed in |33 , by solving the five-dimensional 
field equations in the bulk. The cosmological dynamics on the brane has been analyzed by also including the Weyl 
term, and the relation between the anisotropy on the brane and the Weyl curvature in the bulk has been discussed. In 
these models, it is not possible to achieve geometric anisotropy for a perfect fluid or scalar field - the junction conditions 
require anisotropic stresses on the brane. But in an anti-de Sitter bulk, the solutions can isotropize and approach 
a Friedmann type brane. Bianchi I type brane cosmologies with scalar matter self-interacting through combinations 
of exponential potentials have been studied in '36]. Such models correspond in some cases to inflationary Universes. 
In the brane scenario, as happens in standard four-dimensional general relativity, an increase in the number of fields 
assists inflation. The asymptotic behavior of Bianchi I brane worlds was considered in [23] ■ As a consequence of the 
nonlocal anisotropic stresses induced by the bulk, in the limit in which the mean radius goes to infinity the brane 
does not isotropize and the nonlocal energy does not vanish. The inflation due to the cosmological constant might be 
prevented by the interaction with the bulk. 

The study of anisotropic homogeneous brane world cosmological models l24']-'305 has shown an important difference 
between these models and standard four-dimensional general relativity, namely, that brane Universes are born in an 
isotropic state. For Bianchi typ e I and V geometries this type of behavior has been found both by exactly solving the 
gravitational field equations [23 , or from the qualitative analysis of the model . A general analysis of the anisotropy 
in spatially homogeneous brane world cosmological models has been performed by Coley who has shown that 
the initial singularity is isotropic, and hence the initial conditions problem is naturally solved. Consequently, close to 
the initial singularity, these models do not exhibit Mixmaster or chaotic- like behavior Based on the results of 
the study of homogeneous anisotropic cosmological models Coley conjectured that the isotropic singularity could 
be a general feature of brane cosmologies. This conclusion has been contested as a result of a perturbative analysis of 
the dimensionless shear (38i] . The existence of a decaying mode in scalar perturbations that grows unbounded in the 
past seems to suggest that anisotropy also grows unbounded in the limit of small times [3^. This result is based on 
including, through perturbations, generic inhomogeneities in the cosmological model and these are responsible for the 
unbounded growth of the anisotropy near the singularity. However, in a qualitative numerical study of the asymptotic 
dynamical evolution of spatially inhomogeneous brane-world cosmological models close to the initial singularity, Coley, 
He and Lim ^39] have shown that spatially inhomogeneous G2 brane cosmological models with one spatial degree of 
freedom always have an initial singularity, which is characterized by the fact that spatial derivatives are dynamically 
negligible. From the numerical analysis they have also found that there is an initial isotropic singularity in all of 
these spatially inhomogeneous brane cosmologies, including the physically important cases of radiation and a scalar 
field source. The numerical studies indicate that the singularity is isotropic for all relevant initial conditions. A 
similar result has been obtained by using the covariant and gauge-invariant approach for the analysis of the linear 
perturbations of the isotropic model J^h, which is a past attractor in the phase space of homogeneous Bianchi models 
on the brane p[ol ] . Therefore one can conclude that brane Universes are born with naturally built-in isotropy, contrary 
to standard four-dimensional general relativistic cosmology |4l|. The observed large-scale homogeneity and isotropy 
of the Universe can therefore be explained as a consequence of the initial conditions. 

On the other hand most of the studies regarding the behavior of the anisotropy in brane world cosmological models 
have been done at a qualitative level, and have not provided explicit exact representations for the anisotropy. Also 
many observationally important questions, like, for example, the maximum value of the anisotropy or the moment 
in the evolution of the Universe when this maximum occurred have not been answered yet. It is also not clear for 
what type of equation of state or range of parameters of the cosmological fluid on the brane the initial singularity is 
isotropic, or what is the effect of an inflationary phase on the initial anisotropy. 
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It is the purpose of the present paper to consider the evolution and dynamics of the anisotropy of homogeneous 
anisotropic (Bianchi type) brane world cosmological models in a systematic way. As a first step an evolution equation 
for the anisotropy parameter (describing differences in the time expansion of the Universe along the three principal 
axis) is derived. From mathematical point of view it is a separable first order differential equation for Bianchi types 
I and V, and a Bernoulli type equation for the other Bianchi types. By integrating the evolution equation, one can, 
generally, obtain the anisotropy parameter as a function of the energy density and the volume scale factor of the 
Universe. The study of the behavior of the anisotropy parameter near the singular state shows that generally the 
initial value of this parameter is dependent on the equation of state of the cosmic matter. A high density cosmological 
fluid obeying a barotropic equation of state starts its evolution on the brane from an isotropic geometry, while the 
expansion of a pressureless dust, in an anisotropic space-time, is similar to the standard four-dimensional general 
relativistic one. But for inflationary models the behavior of the anisotropy parameter is similar in both brane world 
cosmological models and standard four-dimensional general relativity. 

The present paper is organized as follows. The basic equations of the brane world cosmological models are presented 
in Section IL The evolution of the anisotropy of Bianchi type I and V models is considered in Section III. In Section 
IV the anisotropy of arbitrary Bianchi type models is analyzed. Finally, in Section V, we discuss and conclude our 
results. 



II. GRAVITATIONAL FIELD EQUATIONS IN THE BRANE WORLD MODEL 

On the 5-dimensional space-time (the bulk), with the negative vacuum energy A5 and brane energy-momentum as 
source of the gravitational fleld, the Einstein field equations are given by 



Gij — k^Tj,/, 



Tij = -A591J + SiY) [-Xgjj + Tf/"-] 



(1) 



In this space-time a brane is a fixed point of the Z2 symmetry. In the following capital Latin indices run in the range 
0, 4, while Greek indices take the values 0, 3. 

Assuming a metric of the form ds^ = {ninj + gjj)dx^dx'^, with nidx^ = dx the unit normal to the x = const, 
hypersurfaces and gij the induced metric on % = const, hypersurfaces, the effective four-dimensional gravitational 
equations on the brane (the Gauss equation), take the form |l5lll6| : 



where S'^i, is the local quadratic energy-momentum correction. 



E„ 



(2) 



(3) 



and Efj^i, is the non-local effect from the free bulk gravitational field, the transmitted projection of the bulk Weyl 
tensor Cjajb, Eij = CiAjsn^n^ , with the property Ejj E^j^i^SjSj as x ^ 0. 

The four-dimensional cosmological constant, A, and the coupling constant, ^4, are given by A = /c| (A5 -I- fc|A^/6) /2 
and fc| — fc|A/6, respectively, with A the vacuum energy on the brane. 

The Einstein equation in the bulk and the Codazzi equation, also imply the conservation of the energy momentum 
tensor of the matter on the brane, Ui/T^" — 0. Moreover, the contracted Bianchi identities on the brane imply that 
the projected Weyl tensor should obey the constraint D^E^J^ = kfD^Sfj," . 

For any matter fields (scalar field, perfect or dissipative fluids, kinetic gases etc.) the general form of the brane 
energy-momentum tensor can be covariantly given as |42| | 



(4) 



The decomposition is irreducible for any chosen 4- velocity . Here p and p are the energy density and isotropic 
pressure, and h^y — g^^^ + u^Un projects orthogonal to u^^ . The energy flux obeys 9^ = g<^i>, and the anisotropic 
stress obeys tt^^ = 7r<^i^>, where angular brackets denote the projected, symmetric and trace-free part: 



(5) 



The symmetry pro per ties of £"^1/ imply that in general we can decompose it irreducibly with respect to a chosen 



4- velocity field as 



I pert! 

m 



E^, = -k^ 



U 



(6) 



4 



where k = ^5/^:4, U is a scalar, a spatial vector and V^i, a spatial, symmetric and trace- free tensor. For 
homogeneous models — and P^j^ = p5| . Hence the only non-zero contribution from the 5-dimensional 
Weyl tensor from the bulk is given by the scalar or "dark radiation" term U. 

In the following we shall consider only homogeneous and anisotropic brane geometries. In particular, those which 
have space-like surfaces of homogeneity (i.e. a G3 acting simply transitively on a V3). The corresponding cosmological 
models fall into nine classes of equivalence, the so-called Bianchi models. It is useful to classify the nine types into 
two disjoint groups, depending on the different properties of the isometry groups (the Lie groups), called class A and 
class B. 

At each point of the space-time on the brane we take a local orthonormal tetrad, where the metric can be written 

as 

where the differential one-forms uj^ will be taken as iv^ = dt and — 0^(^)17% where fi* are time-independent 
differential 1-forms and ai{t), z = 1, 2, 3 are the cosmic scale factors. The time independent 1-forms obey the relations 
dfl^ = — ic^;0'^ A n'', where the c^i are the canonical structure constants and A denotes the exterior product [43l |. 
Although the Ansatz for the metric given by Eq. Q is not the most general one that one can chose, it is sufficient to 
display all the main features of the behavior of Bianchi geometries. 

With the help of the scale factors one can define the following variables: V = n?=i '^i (volume scale factor). 
Hi = ai/a-i^i = 1,2,3 (directional Hubble parameters), H = (l/3)X]i=i^i (mean Hubble parameter) and AHi = 
Hi — H, i — 1, 2, 3. By using the definitions of H and V we immediately obtain H — V/3V, where a dot denotes 
the derivative with respect to the cosmological time t. 

According to the definition of the energy-momentum tensor on the brane, Eq. in the general case of Bianchi 
type geometries, the symmetry of the space-time allows different spatial components of T^^. There are several physical 
processes that could generate an anisotropic energy momentum tensor, with 7^ Tf 7^ Tg , like magnetic fields, heat 
transfer and/or viscous dissipative processes in the cosmological fluid on the brane. The most important of these 
processes are the bulk viscous type dissipative processes, which are the main sources of entropy generation in the early 
Universe. However, the effect of the bulk viscosity of the cosmological fluid on the brane can be considered by adding 
to the usual thermodynamic pressure p the bulk viscous pressure H and formally substituting the pressure terms in the 
energy-momentum tensor by Pef f ~ p + 11 |2H . Therefore the consideration of the bulk viscosity of the cosmological 
fluid does not lead to an anisotropic pressure distribution. The viscous dissipative anisotropic stress tt^j/ of the matter 

on the brane satisfies the evolution equation T2ft.^/iJ^7T^y -I- TTap = —^-qaaf} — rjT {t2U^ /'irfT).^ iTap where is the 

shear viscosity coefficient, T2 = 27^/32, with (32 the thermodynamic coefficient for the tensor dissipative contribution to 
the entropy density and is the shear tensor. Generally, it is assumed that the dissipative contribution from the 
shear viscosity in the early Universe can be neglected, 77 « 0|. Consequently, in the foUowings we consider that 
the anisotropic stresses of the matter on the brane also vanish, « 0. We suppose that in Eq. the heat transfer 
is zero, that is, we take = 0. All these approximations are standard in the analysis of the physics of the very early 
Universe. Therefore in the following we assume that the pressure distribution of the cosmological fluid on the brane 
is isotropic and the fluid pressure satisfies a barotropic equation of state of the form p — p{p). 

For any homogeneous model the conservation equations of the energy density of the matter p on the brane and of 
the dark radiation U can be written as 

p + 3{p + p)H = 0, (8) 



U + iHU = 0, (9) 
leading to a general dependence of p and U of V of the form 

V = ^.U = ^. (10) 



where 

w = exp 

and Co > and Uq > are constants of integration. 



y4/3' 

dp 



p + p{p) 



(11) 
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The modified Einstein gravitational field equations on the brane can be written in the form of the standard Einstein 
four-dimensional field equations, 

G^_iij = klTlfJ-^\ (12) 

where T^a/^^ = -Ag^^/fcl + T^^ + klSf,^ - 18]. Then the effective total energy density, pressure, 

anisotropic stress and energy flux for a perfect fluid are pf"^/^) = A/fcl + p {1 + p/2X) + (p/kfX^U, p^^ff^ = 

p - A/kl + (p/2A)(p + 2p) + {2/klX)U, nl^^^ = {6/kjX)Vf,^ and ql""^^^ = (6fc|A) Q^. Since for homogeneous 
cosmological models = Vfj.u = 0, it follows that in the case of a perfect cosmological fluid there is a close anal- 
ogy between the gravitational field equations on the brane and standard four-dimensional general relativity, with 
the role of the standard energy density and pressure played by p^'^^f'^ and p'^^ff'f , respectively. The formal analogy 
between standard four-dimensional general relativity and brane world cosmology allows the immediate extension of 
the Collins-Hawking definition of isotropization of a cosmological model 44] to the case of brane Universes. 

Hence, we say that a brane world cosmological model approaches isotropy if the following four conditions hold as 
t^oo: i) the Universe is expanding indefinitely and > n) T'^'^ff)'^^ > and t'^'^S m /t'^'^S ^ 0, i = 1,2,3. 
y(e//)Oi^y(e//)oo represents an average velocity of the matter on the brane relative to the surfaces of homogeneity. If 
this does not tend to zero, the Universe would not appear homogeneous or isotropic iii) the anisotropy in the locally 
measured Hubble constant a/H tends to zero, a/H — > and iv) the distortion part of the metric tends to a constant. 
In condition iii) = (l/2)(Tpy(T^'' represents the shear of the normals n^. For a metric of the form ds'^ = dt^ — 
exp {2a) [exp {2j3)]^^ fi^fi", where Q.^ are one-forms that are not exact in general, a is a time dependent function and j3 
is a symmetric traceless matrix, the shear tensor is defined as a^i, — [exp (/3)]^^ [exp (—/?)] + [exp (/3)];^j, [exp (— /S)]^^ 

Q. For Bianchi class A and B models the shear is given by cr^ = a^yO^'' 12 = (1/2) (eLi - ^H'^) ■ 

Spatially homogeneous models can be divided in three classes: those which have less than the escape velocity (i.e., 
those whose rate of expansion is insufficient to prevent them from recollapsing), those which have just the escape 
velocity and those which have more than the escape velocity Q- Models of the third class do not tend, generally, to 
isotropy. In fact the only types which can tend toward isotropy at arbitrarily large times ar e ty pes /, VIIq and 
Vllh ■ For type Vllh there is no nonzero measure set of these models which tends to isotropy Q] . The Bianchi types 
that drive flat and open Universes away from isotropy in the Collins-Hawking sense are those of type VII. 

As an indicator of the degree of anisotropy of a cosmological model one can take the mean anisotropy parameter, 
defined according to 



A 



=1 



For an isotropic cosmological model Hi = II2 — — H and A = Q. The anisotropy parameter is an important 
indicator of the behavior of anisotropic cosmological models, since in standard four-dimensional general relativity it 
is finite even for singular states (for example, A — 2 for Kasner-type geometries 120]). The time evolution of A is a 
good indicator of the dynamics of the anisotropy. 

For a homogeneous brane world model filled with a perfect fluid satisfying a barotropic equation of state the 
conditions i), ii) and iv) of the Collins-Hawking definition of isotropization are naturally satisfied. Hence, if we 
consider only expanding brane worlds, H > and condition i) holds. From the choice of the matter content, and 
due to the symmetries of the energy-momentum tensor, we have t'^'^^^'^^^ = p'-'^if^ > and t'^'^^^^^^ = 0. The choice 
of the geometry implies that condition iv) is also satisfied. Therefore, the value of the parameter a/H is the main 
indicator of the isotropic/anisotropic behavior of a brane world cosmological model. As one can see from Eq. (|13f) 
the quantity a/H is proportional to the square root of the anisotropy parameter, a/H ~ and so, according to the 
Collins-Hawking definition, if A — > 0, a brane world cosmological model will isotropize in the large time limit t — > 00. 

The formal mathematical similarity between standard four-dimensional general relativity and brane world theory 
can be also used to extend the Hawking-Penrose singularity theorem |4^ to brane cosmologies. From the definition 
of the effective energy density p'^'^lf^ and of the effective pressure p^'^fl^ it follows that for a linear barotropic fluid 
with p = (7— l)p, 1 < 7 < 2, the effective energy-momentum tensor on the brane satisfies both the strong and 
weak energy conditions, which can be expressed as {T'^'^ff^t^-' ^ {1/2) gt^^T^'^^^^) Uf,Ut, > and T'^'^f f^f'^u^.u^ > 0, 
respectively, where it^ is an arbitrary timelike four-vector. The first of these conditions implies that the sum of the 
local energy density and pressure is non-negative, p^'^f^^ + p^'^^^'' > and pf"^//) -|- ^p'-'^^^^ > 0. The second condition 
requires that the local energy density be non-negative in every observer's rest frame, p^'^^f ^ > and p^*^/^) -\-p(^ff ) > q. 
Therefore the Bianchi-type brane spacetimes filled with a perfect linear barotropic fluid are singular, since they satisfy 
the following conditions: a) R^i,u^u'^ > for all timelike vectors b) u\^^Ry}^xe[aUp\u^u'^ ^ for a vector tangent 
to some geodesic c) there are no closed time-like curves and d) either i) there is a closed trap ped surface or ii) there 
is a point p for which < for all of the vectors tangent to the past light cone of p |43 . 
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III. EVOLUTION OF THE ANISOTROPY IN BIANCHI TYPE I AND V MODELS 



For a better understanding of the dynamics of the anisotropy in the brane world we consider first in detail the 
evolution of the mean anisotropy parameter A in Bianchi type I and V geometries. From a formal point of view these 
two geometries are described by the line element 



'^^^dz^ 



(14) 



The metric for the Bianchi type I geometry formally corresponds to the case a = 0, while for the Bianchi type V 
case we have a = 1. 

To study the time behavior of the anisotropy parameter A near the initial singular point on the brane we need only 
the Einstein field equations involving the time derivatives of Hi^i = 1,2,3, and which are given by 



1 d . 2a2 kl , 



12^^+3 ^4/3' l'A3. 



(15) 



For a = we obtain the (m), i ^ Q field equations for Bianchi type I geometry, while a — 1 gives the Bianchi type 
V equations on the brane world. 
By summing Eqs. (|15|l we find 



l|(T.i.)^i7 + 3i.-A, ^,/3 
Subtraction of Eq. (|16|l from Eqs. H15|l gives 



2a2 fc| , , fc| 1,4 Uo 



Ki 



i = 1,2,3, 



with Ki, i = 1,2,3 constants of integration satisfying the consistency condition 
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(16) 



(17) 



(18) 



With the use of Eq. l(T7|) the anisotropy parameter defined in Eq. ((T^ becomes 



A 



(19) 



where — (1/3) Y^i=i ^i- Taking the time derivative of Eq. H19fl we obtain the following evolution equation for 
the anisotropy: 



dA H + iH^ , 

— = -2 A, 

dt H 



or, equivalently, 



1 dA 



3i^2 



V [H + 'iH^\=- 



3K^ 



V 



A + + — (p - p) - —PP + -k' . ,., 
Y2/3 '2 ' 12 3 V^4/3 



1.4 ^0 



Taking into account that dp/dV — — {p + p) /V, Eq. (|21|l can be written as 



1 dA _ 2Cl 



1 



A^ dp iK^ {p + p) 
with the general solution given by 

Aip) 



2Cl / A + 2a2Co"^/^«;2/3 + ^ (p _ p) _ i|pp + Ik^UoC'^'^w^/^ [{p + p) w'^]'^ dp - C 
where C is an arbitrary integration constant. 



(20) 



(21) 



(22) 



(23) 
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Eq. H23() gives the general representation of the anisotropy parameter as a function of the energy density of the 
cosmological fluid for the Bianchi type I and V space-times in the brane world scenario. For a cosmological fluid for 
which the thermodynamic pressure p obeys a linear barotropic equation of state of the form p = (7 — 7 = const., 
1 < 7 £ 2, we have p = p^^/V^ , with po ^ a constant of integration, po can be expressed in terms of Co as po — Cq. 
Hence the anisotropy equation Eq. (|23|1 can be immediately integrated to give the general exact dependence of the 
anisotropy parameter on the volume scale factor for Bianchi type I and V geometries: 

A (V) , (24) 

AV2 + 3Q,2y4/3 + kjpoV^'^ + ^pIV^-'^i + k^UnVy^ + C 

where the arbitrary integration constant C 7^ is related, via the field equations, to the constant by the relation 
C = 2K^/3 |23|. The singular state at t = is characterized by the condition V{0) = 0. The value of the anisotropy 
parameter for t = depends on the equation of state of the cosmological fluid. Hence for 1 < 7 < 2, from Eq. H24|) it 
follows 

lim = 0,1 < 7 < 2. (25) 

Therefore the singular state of the high density Bianchi type I and V brane cosmological models is isotropic, with 
A{0) = 0. For the case of the pressureless dust filled anisotropic brane Universes, p = and 7=1. In this case 

The singular state of the dust filled brane Universe is anisotropic, with ^(0) 7^ 0. In the case of the standard 
four-dimensional general relativity (SGR), the behavior of the anisotropy parameter is different from the case of 
brane cosmological models. SGR is recovered if the limits fcs ^ 0, fc ^ and A5 —00 are taken simultaneously. 
Therefore the anisotropy parameter in standard four-dimensional general relativity A^'^^ is given by 

aSGR ^ 21^ (2'7) 

^ ' (Ay2+ 3^2^4/3 +^2^gy2-7 + C7)' ^ ' 



Near the singular state. 



3/<'2 

hm^^««(y) = — >0,V7e [1,2]. (28) 



In fact, one can show that for barotropic matter filled standard four-dimensional general relativistic Bianchi type I 
and V models A'^°^ (0) < 2 jl^. 

Therefore in the case 7 = 1 we obtain the following relation between the initial values of the anisotropy parameters 
A and A^'^^ of the brane world models and of the SGR, respectively: 

^(0) = ^—M- (29) 

In Bianchi type I and V geometries there is also a simple proportionality relation between the shear scalar and 
the anisotropy parameter: 



(30) 



IV. BEHAVIOR OF ANISOTROPY IN ARBITRARY TYPE BIANCHI MODELS 

For arbitrary type Bianchi geometries the (m) ,i 1,2,3 components of the gravitational field equations on the 
brane can be represented in the following general form: 

l|(Fi/.) = (ai, a2, as) + K + ^ {p - p) - ^pp + \k'^. ^ = 1, 2, 3, (31) 
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Bianchi type 


Cl 


C2 


C3 


I 











II 


1 








VIo 


1 


-1 





VIIo 


1 


1 





VIII 


1 


1 


-1 


IX 


1 


1 


1 



TABLE I: Values of the constants Ci, i = 1, 2, 3, for the Bianchi type A models 0- 

where Fi (oi, 02, 03) ,i = 1, 2, 3 are functions which depend on the Bianchi type. For the class A models 

[c2al ~ cgo^)^ - (ciaf)^ 



Fi (ai, 02, as) 



(32) 



where the constants Ci, z = 1, 2, 3 define the Bianchi type, and are given in the table. 

F2 (oi, a2, as) and (ai, a2, as) can be obtained by a cyclic permutation of the elements in the numerator of Fi. 
For Bianchi types V and Ylh 



Fi (ai,a2,a3) 



F2 (01,02,03) 



-2^ + 24^ 



^2"3 



F3 (01,02,03) -2 



of 



(33) 



(34) 



(35) 



with oo, qo and 6 constants. 

If we take qq = b ^ we obtain Bianchi type V, for go : ^ 7^ we obtain Bianchi type VI/i (/i 7^ 0), while for qq = —1 
we have Bianchi type III |4^. 

For Bianchi types IV and Vllft, {h ^ 0) the functions Fi are slightly more complicated. Thus for Bianchi type IV 



Fi (01,02, 03) ^ 



af 2afa2 



(36) 



F2 (01,02, 03) = — 



"3 r2 

2„2 ' n 2-1 ' 
2 ^"2 



2a{a 



(37) 



F3 (01,02,03) 



2a\a 



2„2 

2 



2al- 



with / corresponding to the off-diagonal term [43. For Bianchi type VII^ {h 7^ 0) 

T-l + f^a^ _ oi 
— a\ \ 03 oi 



Fi (ai, 02, 03) = 2A^o^02^ 



02 
02 



F2 (01,02,03) = -Fi (01,02,03) . 



(38) 

(39) 
(40) 



F3 (01,02,03) = 0, 



(41) 
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where A =constant^ 

Generally there is one more field equation, the (00) equation, but it will not be used in the proof of the results. 
By adding Eqs. (|3HI on obtains 

l|(yi/) = H + 3H^=F{a„a„a,)+A+^{p-p)- ^pp + \k'^,, (42) 



where 

3 

3 



1 ^ 

F (01,02,03) ^ ^'^Fi (01,02,03). (43) 



Substraction of Eq. H42|) from Eq. 131|) and integration of the resulting equation gives 

Ag. ^ „. - H ^ + y I- C-) ■ g l'^) . t*-)! ,K . . , 2. 3. (44, 

where 

AF^ (01,02,03) = Fi (01,02,03) - F (01,02,03) ,i = 1,2,3, (45) 

with the property 

3 

^ AF, (01,02, 03) -0, (46) 

i=l 

and Ki, z = 1, 2, 3 are constants of integration satisfying the condition Y^i=i — ^■ 

Therefore for an arbitrary Bianchi type geometry the anisotropy parameter can be represented in the following 
exact form: 



where 

3 

3 



K^ = \j:^Kl (48) 



1=1 



G' 



1 ^ 
27 ^ 



i=l 



H 



-dV 



(49) 



and 



AF,; 

IT 



dV. 



(50) 



Taking the time derivative of Eq. H47I) . and changing the time variable to V ^ it follows that in an arbitrary Bianchi 
type geometry the anisotropy parameter satisfies a Bernoulli type first order differential equation of the form 



dA 
dV 



— ln{K^ + G^ + L) 



A 



2V 



3 (K2 + G2 + L) 
with the general solution given by 

A{V) 



^4 



F (01,02, 03) + A + — (p-p) - -k 



12' 



1 



3 y4/3 



4 ^0 



3 + G2 + L) 



/ ^ (ai, 02, 03) + A + ^ (p - p) - fipp + ifc4Z^oV^"^/3 



(51) 



(52) 
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where C is an arbitrary constant of integration. 

For a brane cosmological fluid obeying a linear barotropic equation of state p = (7 — 1) p, the anisotropy parameter 
for arbitrary Bianchi type cosmological models is given by 

3^27 (if 2 + + i) 



A {V) ^ p , (53) 

/ VF [ai {V) , a2 {V) , 03 (V")] dV + Al/2(i+7) + kjpnV^+'r + %plV^ + /c^^^o 1/2(1/3+7) + cV^-< 

We shall consider in the following the behavior of the anisotropy parameter for arbitrary type Bianchi space-times, 
filled with a linear barotropic cosmological fluid, obeying an equation of state of the form p = (7— l)p, 1 < 7 < 2. For 
this form of cosmological matter the conditions p'^'^^f^ > 0, p'^'^^f^ + p('^ff '> > Q and p'^'^^f^ + Sp^"^^-^' > are satisfied 
on the brane. Therefore the corresponding cosmological models are singular Since in most Bianchi types one 

does not know any exact solution, in order to find the behavior of the anisotropy at early times, it is necessary to 
use some asymptotic solutions obtained, near the singularity, by approximate methods, in the limit of small values of 
the time parameter. In standard four-dimensional general relativity it is concluded that for all Bianchi types there 
exists a Kasner-like "vacuum phase" near the singularity, that is, in general, Einstein's vacuum equations are the 
first order approximation of the equations with a nonvanishing matter term for t — > 0. This idea has been proposed 
a long time ago by Belinskii, Lifshitz and Khalatnikov j^. The general argument comes from the consideration of 
a Kasner type metric ds^ = dt^ — t^^dx^ — t^^dy^ — t^^dz'^, with pi,i = 1,2,3 constants satisfying the conditions 
J2^=iPi = 12^=1 Pi — 1- For a linear barotropic equation of state it follows from the Bianchi identity (jSJ that energy 
density of the matter behaves like p ~ f-7(pi+P2+P3) ^ Therefore for asymptotes fulfilling the Kasner constraints the 
matter term may be neglected in comparison with terms coming from the geometric part of the field equations and 
containing the second time derivatives. A Kasner type behavior near the singularity for a Bianchi type IX brane 
world geometry has been discussed in "2^ 

From Eq. (|53|l one can derive the behavior of the anisotropy parameter near the singular state on the brane for all 
Bianchi types. We suppose that the brane Universe starts its evolution from a singular state, with (0) = 0, i = 1, 2, 3. 
For < > and for an expanding geometry, the scale factors are monotonically increasing functions of time, which for 
small t can be represented in the Kasner form 

a, ^1/"',^ = 1,2,3, (54) 

with nii constants satisfying the conditions < < 1, i = 1, 2, 3 and "t^i — 1. 

The existence of a such a representation for the scale factors near the singular point on the brane has been explicitly 
proven, in the case of Bianchi type I and V geometries, in |27j . The mean Hubble parameter behaves near the singular 
state like H ^V~^. 

With these assumptions it is easy to find the behavior of the integrals involving the functions Fi and Gi. For 
Bianchi class A models, / VF [ai {V) , 02 {V) , 03 {V)] dV can be written, near the singular point, as a sum of terms 
of the form V^"'' , i = 1, 2, 3 and T/2m,+2m, ^i^j- 

3 3 

VF[al{V),a2{V),a3{V)]dV = Y,a^V^"''+ E (55) 

where ai^i = 1,2,3 and f3ij,i,j = 1,2,3 arc constants. 

Thus, for example, for the integral involving the function Fi (ai, a2, 03) we obtain 

/ VFi {ai,a2,a3)dV = ^T/4rn. _ _C2C3 y2{m,+ms) ^ _4_y4m, _ ^yirm _ (gg) 

J 8m2 2 [1712 + ms) 87713 8mi 

For the other Bianchi types, the integral is a sum of terms of the form V'^~^"^\i = 1, 2, 3 or ]/2+2m3-2™i-2m2 g^^^ 
cyclic permutation of this term. Due to the relation mi -I- m2 + — 1 the last expression can be transformed to the 
form V^"^% z = 1, 2, 3. In the limit V" ^ 0, all these terms tend to zero. 

The behavior of the function and L can be analyzed in a similar way, and one can easily show that limy^o G^ = 
limy^o L = Q. Therefore for all Bianchi type cosmological models the anisotropy parameter on the brane has the 
general property limy^o ^ = 0. In the limit of large times, V ^ oo and all Bianchi models (except type VII) 
isotropize, with A — > 0. 



V. DISCUSSIONS AND FINAL REMARKS 



The time evolution of the anisotropy parameter for a linear barotropic fiuid in the brane world model is very different 
from the standard four-dimensional general relativistic case. In homogeneous and anisotropic brane world models the 
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Universe starts from an isotropic state, with A = 0. The anisotropy of the Universe is increasing in time, and reaches 
a maximum value after a finite time ^max- For time intervals so that t > tmax; ^ is a decreasing function of time which 
generally tends, in the large time limit, to zero. In standard four-dimensional general relativity the Universe starts 
its evolution from a singular state with maximum anisotropy and reaches, for all Bianchi models except type VII, an 
isotropic state for i — > oo (44] . For Bianchi type I and V geometries the maximum value of the anisotropy parameter 
is obtained for values of the volume scale factor V satisfying the equation 

2AV + Aa'v'/' + (2 - 7) klpoV'-'' + fpl (1 - 7) V'-'-< + h^U^V-''^ = 0, (57) 

o 

which follows from the condition dA/dV = 0. In the case of the Bianchi type I brane filled with a stiff cosmological 
fluid (7 = 2), and for a negligible small cosmological constant, A = 0, the maximum value of the volume scale factor 
is given by 

The maximum value of the anisotropy parameter is given by 

^max = A (Kiax) - —373 ■ (59) 

klpo + v/f ^ + C 

The value of the cosmological time for which the A reaches its maximum value can be found from Eq. (|15|) , and it 
is given by 

/•Vmax xdX 

The integral in Eq. H60I) cannot be expressed in a simple analytical form. 

If in the early stages of evolution of a stiff (7 = 2) cosmological fluid on a Bianchi type I brane the dark radiation 
term can be neglected, as being negligible small, hit) ~ 0, then the existence of a maximum of A also requires a 
non-zero cosmological constant, A 7^ For Uq = the maximum value of the volume scale factor is given by 

'^s Po 



^-x=(^) , (61) 



and the time necessary for the brane Universe to reach this state is 

xdx 



y^3Aa;4 + Ca;2 + MPo 
The maximum value of the anisotropy parameter is, in this case. 




(62) 



^max = A (Knax) = 7 ' (^3) 

C 



PO {kl + ^kl) 



This type of behavior of A, with an initial monotonically increasing evolution from zero up to a maximum value, 
followed by a decrease to zero, specific to brane world cosmological models, is a direct consequence of the presence, 
in the gravitational field equations, of the term quadratic in the energy density. 

In the present paper we have considered in a systematic manner the time evolution of the anisotropy parameter A 
in the framework of homogeneous, arbitrary Bianchi type brane world cosmological models, with the matter content 
consisting of a perfect barotropic cosmological fluid. For Bianchi type I and V exact representations of A can be 
obtained from the field equations and therefore the initial behavior of A can be explicitly derived from the study of 
the exact representation near the singularity. In order to obtain the form of A at the initial moment of the cosmological 
evolution for the other Bianchi types we have used the crucial assumption that near the singularity the metric of the 
brane world is of Kasner type. For a "normal" matter filled brane world, satisfying a linear barotropic equation of 
state, the behavior of A is very different from the standard four-dimensional general relativistic case. In this case 
the anisotropy parameter has the remarkable property limt_+o A{t) = liuit^oo -^{1) = 0, in sharp contrast to the SGR 
case. However, in the case of pressureless dust, the initial values of the anisotropy parameter are identical in both 
brane world and standard four-dimensional general relativistic cosmological models. 
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